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ABSTRACT 

One of the possible origins of short gamma-ray bursts (SGRBs) is merging of compact binaries, and the effect 
of large kick velocity is a signature that can be used as an observational test for this hypothesis. Intracluster 
SGRBs that escaped from a host galaxy in a galaxy cluster are interesting in this context, since they would 
escape more easily by cluster tidal force, and would have brighter afterglow luminosity by dense intracluster 
gas, than those in general field galaxies. Here we calculate the escape fraction of compact binaries from their 
host galaxies in a galaxy cluster, and discuss some observational implications. We found that the escape fraction 
strongly depends on the nature of dark matter subhalos associated with member galaxies. If the amount of dark 
matter around member galaxies is not large and the gravitational potential for an escaping binary is determined 
mostly by stellar mass, most of SGRBs should escape and be observed as hostless, which is a much higher 
fraction than those in the field. Hence, statistics of intracluster SGRBs could give important information about 
the dark matter distribution in galaxy clusters, as well as hints for the origin of SGRBs. 

Subject headings: binaries: general — gamma rays: bursts — galaxies: clusters: general 
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1. INTRODUCTION 

Short gamma-ray bursts (SGRBs) occur not only in star 
forming galaxies but also in early type galaxies, indicat- 
ing a longer time sc ale of the delay from star formation 
to the SGRB events fcehrels et alH 120051: IVillasenor et al. I 
120051 iBerger et al. 1 120051) than that of the long gamma-ray 
burs ts (LGRBs) th a t are ass o ciated with massive star collapses 
rseelPiranl d2004l) : iNakarl (120071) and Lee & Ramirez-Ruiz I 
for reviews]. A strong candidate for the origin of 
SGRBs satisfying this nature is merging of compact object 
binaries (e.g., double neutron star binaries or neutron star - 
black hole binaries). 

Compact binaries are expected to have large systemic ve- 
locities by the kick given at the time of compact object forma- 
tion, as inferred from the large proper motion of pulsars. Be- 
cause of the long time scale and the large velocity, SGRBs are 
expected to occur with a large offset from the centers of galax- 
ies compared with stellar distribution. There are some model 
calculations of this effect (Bulik, Belczyhski, & Zbijewski 
1999; Bloom, Sigurdsson, & Pols 1999; Fryer, Woosley,& 
Hartmann 1999) and such an expectation is indeed consis- 
tent with observations, though the statistics is still limited 
(Lee et al. 112001 iFox et al. ll20lEllBloom"et1a]l2007h . How- 
ever, truly intergalactic SGRBs that are ejected from their host 
galaxies would be difficult to identify, because of the difficulty 
to identify their host galaxies. Furthermore, we do not expect 
a strong afterglow emission in low density environment in the 
intergalactic field, making the id e ntification eve n more diffi- 
cult [see also lBloom etaiTd2006h : lNakar I d2007l) l. 

SGRBs in galaxy clusters are particularly interesting in this 
context, since compact binaries would be more efficiently 
stripped away from their host galaxies and become intraclus- 
ter SGRBs, because of the strong tidal force in the cluster 
gravity field. The idea of tidal st ripping is consistent with ob- 
servations of in t raclus t er stars (IVflchez-Gomez et al. Ill994t 
Okam ura etaTI 120021; [D urrell et al. I l2002t iGal-Yam et al] 
20031; iGerhard et al. I l2007t) . We expect that such intraclus- 



1 Department of Astronomy, School of Science, Kyoto University, Sakyo- 
ku, Kyoto 606-8502, Japan 

2 niinou@kusastro.kyoto-u.ac.jp 



ter SGRBs could have brighter afterglows because of the 
denser environment of the intracluster gas than that in the in- 
tergalactic field, which would help the accurate determina- 
tion of SGRB locations. A few SGRBs have already been 
observed in galaxy clusters, though the sa mple is still small 
( IBerger et al. 1120071: IShin & Berger1l2007h . The statistics of 
such intracluster SGRBs in the future data would then give us 
some information about the origin of SGRBs, and the gravi- 
tational potential of galaxies in a cluster. 

Here we calculate the expected event rate of such intraclus- 
ter SGRBs, assuming that SGRBs are indeed the products of 
compact binary mergers. We also discuss possible implica- 
tions for future observations. We assume ACDM cosmology 
with fl A = 0.7, il M = 0.3, and H = 70 km s" 1 Mpc" 1 . 

2. FORMULATIONS 

2.1. Models of a Galaxy Cluster and Its Member Galaxies at 

z = 

First we construct a model of a galaxy cluster at z = 0, 
and will consider its evolution in the next subsection. As a 
typical galaxy cluster, we consider a spherically symmetric 
cluster with a total cluster mass of 1O 15 M , a virial radius of 
R wil . = 1 Mpc, and the Navarro-Frenk- White (NFW, Navarro, 
Frenk, & White 1995) density profile for dark matter (DM) 
with a concentration parameter of Chaio = 5. All these pa- 
ramete rs are typical for a rich cluster [e.g.. Sch midt & Allen I 
(2006)]. We assume the Schechter function for the luminosity 
function of member galaxies, with the shape parameters (i.e., 
L* and a) independent of radius from the cluster center, R. We 
use the estimate of the Schechter paramet ers by the Sloan Dig- 
ital Sky Survey data for galaxy clusters (Ha nsen et al. 1120051 
Table 2). The ^-dependence of the normalization of the lu- 
minosity function, 0*, is assumed to be the NFW profile, but 
with the concentration param eter c ga i aX v = 1 -2 that is different 
from that for the dark matter (Han sen et al. 112005b . 

We then determine the stellar mass and density profile 
within a member galaxy as follows. It is well known that 
early-type galaxies are the dominant population in ri ch galaxy 
clusters (lOemler I fl974l iThomas & Katgerdll2006l) . and we 
assume all cluster galaxies to be elliptical galaxies for the sim- 
plicity. The observed luminosity function is in the r' band, and 
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we translate it into the B band by a typical color of elliptical 
galaxies, B— r 1 = — 1 .4. We assume the ratio of the stellar mass 
to light to be 9.8Mp l /L Pl in the B band. These q uantities are 
from the model used by Totani & Yoshii (2000), which can 
reproduce the observed galaxy properties well. Stellar density 
profile withi n member galaxi es is assumed to be the Hern- 
quist profile (Hernquist 1990), with the characteristic radius 
parameter chosen so that a half of the total mass is included 
within the observed half-light radius, rhaif, which is calculated 
by th e power-law relation to B luminosity fitted to observa- 
tions dTotani & Yoshii 1120001) . 

On the other hand, the density profile of dark matter as- 
sociated with the cluster member galaxies as substructure in 
the whole cluster is uncertain and poorly known. Especially, 
the dark matter subhalos extending to the outer region of 
member galaxies are expected to be vulnerable to the tidal 
forces by galaxy interactions and/or the overall cluster gravi- 
tational potential. Here we consider the two extreme cases: (i) 
there is no significant DM substructure or subhalos associated 
with member galaxies, and the gravitational potential well of 
galaxies is determined simply by the stellar mass profile (the 
no subhalo case), and (ii) the DM subhalos are associated with 
member galaxies with a similar amount to the field galaxies 
(the preserved subhalo case). The amount and nature of sub- 
halos depend on the cluster forming processes, and probably 
the reality is between these two extreme cases. 

For the preserved subhalo case, we calculate the virial mass 
(M v j r ) of the subhalos from the stellar mass of a galaxy and 
the universal ratio of the dark-to-baryonic matter Q,M/^b = 
5.9 (ISpergel et al. 1 12003). The virial radius (r vir ) of the 
subhalo is calculated from M v i, and the 1 -dimensional cen- 
tral velocity dispersion (cr v ) of the galaxy, as GM v i r /2r v i r = 
3a„. Here, the velocity dispersion is that for stars calculated 
from the galaxy luminosity us ing the Faber-Jackson relation 
(Ide Vau couleurs & Olson Ifl982l) . An almost similar velocity 
dispersion is obtained also from the stellar mass and rhaif cal- 
culated above and assuming the virial relation, giving a con- 
sistency check for our treatment. We assume the NFW profile 
for the DM subhalos, and the c oncentration paramet er is cal- 
culated by the formula given by lBullock etaLl (1200 ll) for sub- 
halos included in larger virialized halos, which are based on 
cosmological N-body simulations: c ~ c»(M v i r /M* ) r at z = 
where c* = 7, T = -0.3, andM* ~ 1.5 x \Q n lT l M & . 

2.2. Cluster and Galaxy Evolution 

It is well known that most stars in elliptical galaxies formed 
at high redshift (z > 1) and they evolved passively to the 
present time (e.g., Yamada et al. 2005). Major mergers could 
change drastically the stellar mass distribution and gravita- 
tional potential of member ga l axies, but a recent numerical 
simulations by iMurante et al. I (120071) indicates that the ma- 
jority of member galaxies in a cluster do not undergo ma- 
jor mergers, except for the brightest central galaxy which is 
formed in the collision of many galaxies. Therefore we make 
a reasonable assumption that the stellar mass distribution does 
not evolve in cluster member galaxies. 

Though galaxies form at the early epoch of z > 1, the estab- 
lishment of the overall cluster potential should be significantly 
later according to the standard picture of hierarchical structure 
formation. We estimate the epoch of cluster f ormation using 
the ex tended Press-Schechter approximation (lLacey & Cole I 
119931) . It predicts that about half of mass of a 1O 15 M clus- 
ter at z = is already included in the largest progenitor at 
Zf ~ 0.45. Therefore we assume that the tidal force by the 



cluster potential starts to affect member galaxies at this red- 
shift. We do not take into account the evolution of the clus- 
ter potential at z < Zf, and this is a reasonable approximation 
because the time scale of cluster potential evolution is much 
larger than the orbital period of compact binaries in member 
galaxies. The evolution of DM subhalos is difficult to predict 
without detailed numerical simulations, and we simply ap- 
ply the above two extreme models with no evolution, which 
would cover the realistic evolution. 

2.3. Escape of Compact Object Binaries 

In a galaxy with a given luminosity, we can calculate the or- 
bit of compa ct ob ject binaries in the gravitational potential as 
modeled in $2. II if the initial location and velocity are given. 
We calculate the initial velocity by the sum of the original 
stellar velocity v 01g at the location and the kick velocity Vkick 
given when the compact objects are formed. The Maxwell 
distribution having the 1 -dimensional velocity dispersion <r v 
defined in the previous section is assumed for v 01g . The di- 
rection of both v org and Vkick are assumed to be isotropic and 
random. 

Although the velocity distribu tion of observed s ingle pul- 
sars is well fitted by a Gaussian (Hobbs et al.n i2005h . that for 
compact object binaries is not well known. A few observed 
binar y pulsars have bulk motion velocities of ~ 100-200 
km/s dWex, Kalogera. & Kramer l|2000t iRansom et al. 112004 
I Willems. Kaloeera & Henninger ||2004). iBulik et al.l (1 19991) 
Bloom et al.ldl999l) and lFrver et al. I(ll999l) theoretically esti- 
mated the velocities of compact binaries that remain gravita- 
tionally bound after supernova explosions, to be several hun- 
dreds km/s. We assume the distribution of Vkick to be a sin- 
gle isotropic Gaussian, i.e., each 1-dimensional component of 
Vkick is a Gaussian with the standard deviation <t>. We cal- 
culate the cases of two different values of = 100 and 300 
km/s. The mean (the standard deviation) of the corresponding 
I Vkick | distribution then becomes 160 (66) and 480 (200) km/s 
for <j\ =100 and 300 km/s, respectively. 

We solve the motion of binaries until their merger time (the 
time from a compact binary formation to its merger by grav- 
itational wave radiation). Though the merger time generally 
ranges more than three orders of mag nitude, 10 7 " 10 yr (e.g., 
Tutukov & Yungelson 1994; Bulik et al. 1999), we are inter- 
ested in SGRBs in galaxy clusters. Most of galaxies in clus- 
ters have formed their stars at high redshift (z > 2), and ob- 
served SGRBs are typically at z ~ 0.2. Therefore, the time be- 
tween these epochs, i.e., ~ 10 10 is appropriate for the merger 
time in this work. 

We consider that a binary has escaped from its host galaxy 
once its distance from the host galaxy center becomes larger 
than the tidal radius r t id e of the host galaxy after the forma- 
tion of the galaxy cluster, i.e., z < Zf- Here, the tidal radius 
is defined as the radius where the tidal force by the overall 
cluster potential is the same as the binding force in the host 
galaxy. The tidal radius is numerically calculated for a given 
set of galaxy luminosity and location in the cluster. We then 
calculate the mean escape fraction as a function of R, taking 
a weighted average over the host galaxy luminosity, initial lo- 
cation in the host, and kick velocity. 

It should be noted that some fraction of stars are distributed 
at r > r t ;de with the assumed stellar mass profile in host galax- 
ies. Such stars would be stripped from host galaxies and be- 
come intracluster stars. Compact object binaries in such stel- 
lar populations would all contribute to the intracluster SGRBs. 
We find that this fraction is about 2 % in the preserved sub- 
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halo case and 10 % in the no subhalo case. It seems that the 
no subhalo case is preferred (see 14.2b . from a comparison of 
these values with the various observational estimates of the 
abundance of intracluster stars in galaxy clusters. 

3. RESULTS 

The results are shown in Figure Q] where the mean es- 
cape fraction within a given radius from the cluster center, 
/esc(< R), is plotted. We find that the escape fraction largely 
depends on the existence of the DM substructure; the escape 
fraction is modest with / esc ~ 0.2 in the preserved subhalo 
case, while most binaries will be ejected in the no subhalo 
case. The dependence on the radius from the cluster center 
or on the kick velocity is not as significant as the effect of 
subhalos within the parameter ranges investigated. 

For comparison, we calculate the case of field galaxies not 
in clusters. Again, we only consider elliptical galaxies, and 
their proper ties are calcula t ed in the same way for a given 
luminosity. Bullock ^t al. I d2001l) found that isolated halos 
have a different relation between the virial mass and the con- 
centration parameter from that for subhalos, and we adopt 
e* = 9 and T = -0.13 here based on their results. The escape 
fraction averaged over galaxy luminosity is simply calculated 
without taking into account the tidal force by external gravity 
fie ld. Here we use the l uminosity function shape parameters 
of iBlanton et al. I (1200 ll) for field galaxies; though these are 
derived for all types of field galaxies, the luminosity function 
for each galaxy type is rather uncertain. The results are shown 
in Table 1, and we found that the escape fraction in the field 
is not much different from that in galaxy clusters in the pre- 
served subhalo case, while a large enhancement of the cluster 
/esc is predicted in the case of no subhalos. 

4. DISCUSSION 
4.1. Detectability of SGRB Afterglows in Clusters 

Detection of afterglows is necessary to locate a GRB ac- 
curately enough with respect to a host galaxy. We discuss 
here the detectability of a typical SGRB afterglow in the in- 
tracluster medium following th e standard afterglow model of 
ISari. Piran. & Narayan I {1998). We simply use this isotropic 
model without jet structure by using the isotropic equivalent 
total energy. The jet break may reduce the expected flux at a 
later time compared with the calculation here, but this crude 
estimate is sufficient here for our purpose. 

The isotropic-equivalent total energy in gamma-rays of 
SGRBs is distributed in a wide range of E 1 ^ ~ 10 49 - 
10 51 erg, and the total initial kinetic ener gy of the exter- 
nal shock (-Ejso) is ex pected to be simi lar (|Fox et al. 1 120051; 
Soder berg et al. 1120061) . According to Pana itescu & Kumar] 
(2001), we adopt the following parameters: the fraction of 
energy density in magnetic field e# = 10~ 2 4 and that in non- 
thermal electrons e e = 10~ 1,2 ; although these are for LGRBs, 
the values inferred for available SGRBs are not much differ- 
ent. We also assume the power index of the luminosity decay, 
a = -1 where F v oc t a . The typical particle density of intra- 
cluster medium is n ~ 10~ 3 cm" 3 within a few hundreds kpc 
from the cluster center (e.g., Lewis, Buote, & Stocke 2003). 
We assume a typical distance for SGRBs, z = 0.2. 

Then the model predicts the expected flux at the observed 

frequency of v as F v ~ 1.5 nj^ ^^^(f/lO 3 s)" 1 /iJy, 
where £50 = Ei so /(lQ 5Q erg), n_3 = n/(10~ 3 cm" 3 ) and ^15 = 
v/(l0 15 Hz). In the X-ray band (1 keV), the typical flux is 
then vF v ~ 9.3 x 10" 14 erg cm" 2 s _1 at t = 10 3 s, which can 



be detected by existing X-ray satellites (e.g., Gehrels et al. 

2005) . Swift XRT can locate afterglows with an accuracy of 
a few arcsec, and this is reasonably accurate to discuss the 
association of an afterglow with a galaxy at z ~ 0.2. In the 
optical (R) band, this flux corresponds to ~26 mag (AB) at t = 
10 4 s. Most afterglows with this level of brightness have been 
missed in the past and current GRB follow up observations. 
However, afterglows of brightest SGRBs {E lso ~ 10 ) may 
be detectable. 

On the other hand, SGRBs ejected far from their host galax- 
ies in the normal field would occur in much lower-density 
environment. The typical density in general intergalactic 
medium would be estimated as n ~ p c ili,/m p ~ 10~ 7 cm -3 , 
where p c and m p are the critical density of the universe and 
the proton mass, respectively. The expected afterglow flux of 
an intergalactic SGRB is then more than one order of magni- 
tude fainter than those of the faint est SGRB afterglows ever 
observed, suc h as GRB 050509 B dGehrels et al. H2005I) and 
GRB 05091 1 dPage et al. 1120061 upper limit only). Therefore 
it seems difficult to detect an afterglow of such an event. 

4.2. Comparison with Observations 

iBerger et ail d2007l) examined all 16 SGRBs that were fol- 
lowed up by X-ray observations with XRT of Swift or Chan- 
dra, and found that three SGRBs are likely to be associated 
with galaxy clusters, suggesting that the fraction of SGRBs 
in galaxy clusters is about 20 %. Considering the statistical 
uncertainty, it is consistent with the fraction of all stellar mass 
in the universe bound in galaxy clusters (~ 10 %, Fukugita, 
Hogen, & Peebles 1998). 

Among the three SGRBs in galaxy clusters discussed in 
iBerger et all d2007l) . GRB 050509B is apparently associ- 
ated with the likely host galaxy at the cluster redshift of 
Z = 0.226. The offset of GRB 050509B from its host is 21- 
56 kpc, corresponding to 6-16rh a if- The galaxy cluster that 
contains the host is composed of two subclusters, and the host 
galaxy is located at the center of the minor subcluster that 
is a bout 270 kpc away from the center of th e major subclus- 
ter dGehrels et al. Il2005t iBloom et alj|2006l) . The location of 
GRB 05091 1 is in a cluster of z = 0. 165, but its aft erglow was 
too fa int to associate it with any particular galaxy (Pag e et aL~l 

2006) . GRB 050813 has three candidate host galaxies near its 
Swift XRT location, and the galaxies b elong to two differen t 
galaxy clusters at z = 0.72 and z = 1 .8 dBerger 1120051 120061) . 
Clearly, the current sample is too small to derive any implica- 
tions from a comparison with our results. Future satellites for 
GRB study might detect SGRBs more efficiently leading to a 
much larger sample of SGRBs. 

Observations of intracluster diffuse light, stars, type la 
supernovae, and planetary nebulae indicate that some stars 
in a galaxy cluster are in intra cluster medium, perhaps re- 
moved from member galaxies (Vflchez-Go mez et al. 1 1 19941: 
Okam ura et all 120021; iDurrell et al. I 120021: iGal-Yam et alj 
20031 iGerhard et al. 1120071) 7 The fraction of intracluster stars 
in all stars in a cluster is uncertain, but observational estimates 
are typically ~ 5-20 %. These fractions are consistent with 
our estimate of stars outside the tidal radius (see jj2.31 l. On the 
other hand, / esc of SGRBs could be much higher than these, 
up to ~ 80 % depending on the model parameters. If such a 
higher fraction of intracluster SGRBs than that of intracluster 
stars is observed in the future, it would indicate the effect of 
kick velocities on compact binaries, giving a further support 
for the compact binary hypothesis of SGRBs. Note that this 
test is difficult in the intergalactic field, since we do not know 
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the fraction of intergalactic stars and afterglows of intergalac- 
tic SGRBs are difficult to detect. 

5. CONCLUSIONS 

We investigated the escape of compact binaries from their 
host galaxies in galaxy clusters, which is enhanced by the 
tidal force of the cluster gravity compared with general fields. 
We found that the escape probability heavily depends on the 
uncertain distribution of subhalos associated with member 
galaxies. If the DM substructure has been destroyed by in- 
teractions in a galaxy cluster and the escape of a binary is de- 
termined mainly by gravity of stellar mass, most of compact 
binaries in galaxy clusters should escape and become hostless 
intracluster SGRBs. On the other hand, if the DM subha- 
los are associated to member galaxies with a similar amount 
to field galaxies, the enhancement of escape fraction is only 



modest compared with field galaxies: about 20 % for clusters 
while ~ 10 % for field galaxies. 

Though the current observed data set is not sufficient to be 
compared quantitatively with our results, statistics of SGRB 
association with cluster galaxies in the future data will give 
us important information for the dark matter distribution in 
clusters, intracluster stars in clusters, as well as the origin of 
SGRBs. 
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TABLE 1 

ESCAPE FRACTIONS IN GALAXY CLUSTER /esc(< i?vir) AND IN FIELD 



o> (km/s) 


Cluster (i) 


Cluster (ii) 


Field 


300 


0.79 


0.20 


0.10 


100 


0.69 


0.16 


0.09 



NOTE. — The colum ns labeled as Cluster (i) and (ii) give the escape fractions corresponding to (i) the no DM subhalo case and (ii) the preserved DM subhalo 
case discussed in j|2.1| 
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presumed DM subhalo, 0^=300 km/s 

— - presumed DM subhalo, 0^=100 km/s 
no DM subhalo, 0^=300 km/s 

— - no DM subhalo, 0^=100 km/s 



0.5 



500 1000 

R [kpc] 

FIG . 1 . — The escape fraction of compact binaries from host galaxies in a galaxy cluster averaged within R, where R is the distance from the center of the galaxy 
cluster. (The virial radius of the cluster is 1 Mpc.) The models with no dark matter subhalos is shown by thin lines, while the models with preserved subhalos 
are shown by thick lines. The solid and dashed lines are for different values of the standard deviation of the kick velocity distribution, cr^ = 300 and 100 km/s, 
respectively. 



